Asymptotic normalization coefficients and continuum coupling in mirror nuclei 
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Background: An asymptotic normalization coefficient (ANC) characterizes the asymptotic form of a one-nucleon 
overlap integral required for description of nucleon-removal reactions. 

Purpose: We investigate the impact of the particle continuum on proton and neutron ANCs for mirror systems 
from p- and sd-shell regions. 

Method: We use the real-energy and complex-energy continuum shell model approaches. 

Results: We studied the general structure of the single-particle ANCs as a function of the binding energy and 
orbital angular momentum. We computed ANCs in mirror nuclei for different physical situations, including 
capture reactions to weakly-bound and unbound states. 

Conclusions: We demonstrated that the single-particle ANCs exhibit generic behavior that is different for 
charged and neutral particles. We verified the previously proposed relation [l| between proton and neutron 
mirror ANCs. We find minor modifications if the spectroscopic strength is either localized in a single state or 
broadly distributed. For cases when several states couple strongly to the decay channel, these modifications may 
reach 30%. 

PACS numbers: 21.10.Jx, 21.10.Sf, 21.60.Cs, 24.10.Cn 
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I. INTRODUCTION 

The ANC method [2( has proved useful as an indirect 
tool to determine direct capture reaction rates 0-Q both 
to well-bound and weakly-bound states. The ANC of the 
virtual proton decay of a nucleus a -> b + p is related 
to the astrophysical 5-factor for a proton capture reac- 
tion b(p,^f)a at stellar energies, and can be obtained in 
transfer reactions that offer higher cross-sections; hence, 
it provides an alternative to direct reactions at relevant 
energies in astrophysical processes 0, Oj @ • The method 
depends on the peripheral nature of low-energy capture 
reactions, in which case the cross-section is determined 
by the tail of the radial overlap integral between the wave 
functions of the final nucleus and initial colliding systems. 

The ANC, or the nuclear vertex constant, for the de- 
cay a — > b + c is proportional to the asymptotic behavior 
of the wave function representing the relative motion of 
particles b and c. This quantity is closely related to the 
reaction amplitude [5[ and - unlike the spectroscopic fac- 
tor - is invariant under finite-range unitary transforma- 
tions of the nucleon-nucleon (NN) interaction [7[ ; hence 
is less dependent on the choice of a short-range potential. 

Recently, it has been advocated [l| that the charge 
symmetry of the nuclear force could be used to relate the 
ANC of the proton decay to the virtual neutron decay of 
the mirror nucleus. This observation opens an attractive 
possibility to learn about the decay width of hardly ac- 
cessible states in proton-rich nuclei from transfer reaction 



studies in mirror-bound systems using stable beams. It 
remains, however, an important question to what extent 
the charge symmetry argument is sufficient to extract re- 
liable information about the width of an unbound state 
of some nucleus from the ANC in its bound mirror part- 
ner. It is well-known, for example, that threshold effects 
lead to striking differences in the energy spectra of mir- 
ror nuclei having different particle emission thresholds 
Q. Indeed, for near-threshold states, the configuration 
mixing involving scattering states strongly depends on 
the positions of particle emission thresholds in mirror 
systems (the binding energy effect) Q, and on different 
asymptotic behavior of neutron and proton wave func- 
tions. The latter effect leads to the universal behavior of 
cross sections [10|, [H| and overlap integrals fl2l - [l4J in the 
vicinity of the reaction threshold. 

The main objective of this study is to verify the conjec- 
ture of Refs. [lEl-GII that the ratio K of ANCs for mir- 
ror pairs is both approach-independent and interaction- 
independent. To investigate the proposed link between 
proton and neutron mirror ANCs [l], we employ the 
framework of the nuclear Shell Model (SM) for open 
quantum systems (OQS), i.e., the continuum shell model 
(CSM) Qjl, [2fJ, which offers a realistic treatment of the 
configuration mixing in well-bound, weakly-bound, and 
unbound states. 

The paper is organized as follows. In Sec.|TI]wc discuss 
basic features of ANCs. In particular, we review proper- 
ties of ANCs for proton wave functions in bound nuclei 



and the limiting behavior of ANCs for charged particles 
and neutrons. The CSM results for the ANC in mirror 
p- and (sd)-shell nuclei are discussed in Sec. IIIII We em- 
ploy the complex-energy Gamow Shell Model (GSM) and 
the real-energy Shell Model Embedded in the Continuum 
(SMEC) to describe ANCs for bound states and reso- 
nances. We also analyze the dependence of the ANCs in 
mirror systems, and their ratios, on the strength of the 
continuum coupling and configuration mixing. Finally, 
Sec. IIVI summarizes the results of our work. 



For r 3> R, where R is the nuclear radius, ugj is given by 
its asymptotic form 



u ij( r ) = /%W-fj,M-i/2(2«T), 



(5) 



where fyj = j3ij{S c ) is the single-particle ANC 
(SPANC). Therefore, far from the region of nuclear in- 
teraction, I§ c .gj behaves as: 



,1/2 



W) ~ z s 4 fajW^ i+1/2 (2 Kr ). 



(6) 



II. BASIC FEATURES OF THE ASYMPTOTIC 
NORMALIZATION COEFFICIENT 

Let us consider the radiative capture reaction b + c — > 
a + 7, and define the radial overlap function Ig for a 
process a — > b + c. In the asymptotic region, Ig c can be 
written as: 



4W) ~ ;C-«W-^+i/a(2«r), 



(1) 



where W- n j+i/2 is the Whittaker function and Cij is 
the ANC, a quantity characterizing the virtual decay of 
a nucleus into two particles b and c. In Eq. (|Tj). r is the 

relative distance between b and c, k = y2fiS c a /K 2 where 

Sc is the separation energy of particle c in the nucleus 
a, and rj = Zi,Z ' c e 2 '[i / 'ft? 'k where /i is the reduced mass 
of b + c. The quantum numbers £ and j are the orbital 
angular momentum and the channel angular momentum, 
respectively. 

In this paper we shall assume that the particle c is a nu- 
cleon (proton or neutron); hence, Sn,p is the one-nucleon 
separation energy The corresponding radial overlap in- 
tegral can be written as: 



hc-ij\ r ) — 



v J-^ Yl^ a \\4 J (BW J b b )(rmuB 



)(2) 



where ajAB) is a creation operator associated with the 
single-particle (s.p.) basis state |ug) and {r£j\us) is the 
radial s.p. wave function. The sum in ([2]) runs over the 
complete s.p. basis. The squared norm of the radial over- 
lap integral $Z$ defines the spectroscopic factor Sij. 

In the general case of multi-channel coupling, the ANC 
for a radiative capture reaction is defined in terms of the 
Hermitian norm \C\ of all the contributions correspond- 
ing to different couplings of the target state and the state 
in a parent nucleus: 



\C\ 




(3) 



For bound states, Iu c .(a can be well approximated by 



the product of the spectroscopic amplitude SA and the 



1/2 



s.p. radial wave function U£j/r at a s.p. energy —S, 



(a). 



Hence, /3gj is directly related to the ANC [211 ] : 



C, 



u 



s. 



1/2, 



(7) 



The relations (|4|- also hold for many- body resonances. 
Indeed, as demonstrated in GSM studies [T^], the overlap 
function I^ c . ( a is well approximated by the product of the 

1/2 

spectroscopic amplitude SA and the s.p. resonance wave 
function of the average potential, which reproduces the 
Q-value of the reaction studied. Let us also recall that 
the astrophysical 5(, c -factor, in the limit of zero center- 
of-mass energy, E cm = 0, is simply proportional to /3|- 



3] 



5 6c (0)~/3f,. =P%{ri). 



(8) 



A. General properties of SPANCs for charged 
particles 

Let Rf (Rf 3> R) be the radius of the external region 
where the nuclear part of the potential is practically zero. 
In this region (r > Rf) the s.p. wave function ugj is given 
by Eq. ([5]), i.e., it is given by W- v j + i/2(2Kr). For the 
normalized bound state, this implies that A/int+A4xt = 1, 
where A/i n t is the norm of the internal part of the wave 
function, and 






A4xt=4/_ \W^ J+1/2 {2Kr)\ 2 dr (9) 



is the norm of the external part. Provided that the en- 
ergy dependence of A/l n t is weak (which is a reasonable 
assumption even if separation energy is close to zero), 
Pij should strongly depend on the value of the integral 
in Eq. ©. 

In terms of the outgoing Coulomb wave function 
H + (£,r]k,kr), the Whittaker function can be written as 

W_ Vti+1/2 (2Kr) = H+(£ lVkl kr)e^ + ^-^^\ (10) 

where k = in, rj k = —ir], and a(£, rf) is the Coulomb 
phase shift. To discuss the limiting cases, is useful to 
introduce the complex turning point 



Zt 



(v+ vV + ^ + 1)) (ii) 



(4) at which d 2 H+{£, r ]k ,z)/dz' 2 



0. 



According to the standard properties of the Coulomb wave functions [2^,|23J|, one obtains the asymptotic expressions: 

1/4 



it: 



1/2 (2kx) 



r, lt -\ exp [rj - r]\n(rj) - 2\j2r\K,r) for nr <C \z t \ , r] -> +00, 

W_ ri £ + i/ 2 ('2Kr) ~ exp (— kj — J7ln(2«r)) for Kr 3> |z t |. 



(12) 
(13) 
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Near-threshold limit of SPANCs for charged particles 



Let us first consider the limit (fl2j) of 77 — >• +00, which 
corresponds to very small separation energies. Since 77 = 
Ko/k, where kq = ZbZ c e 2 fj,/H 2 , the asymptotic part of the 
s.p. wave function near the particle emission threshold is 



Uij(r) ~/% 



KqT 

2'/ 



1/4 



^3) exp (77 — ?7ln(?7) — 2v2kot) 



(14) 

As the external norm (|9|) must be finite, and exp [77 — 
r)]n.(rf)]/r) 1 ' 2 — > for 77 — > +00, in the limit of very weak 
binding fiij must exhibit the universal 77— dependence: 

ftifa) = WjV /2 exp(7 ? ln77-77), (15) 

where Np(£,j) is a prefactor that depends on the struc- 
ture of the s.p. state, in particular £. To assess how 
quickly the limit (|15jl is reached, we performed calcula- 
tions for the lsi/2, Opi/2, and Od^/2 single-proton states 
in 17 F. The s.p. radial wave function ugj was calcu- 
lated using the Woods-Saxon (WS) potential with the 
following parameters: the strength of the spin-orbit term 
Vso = 3.68 MeV, radius i?o = 3.214fm, and diffuseness 
d = 0.58 fm. We took the Coulomb potential of a spher- 
ical uniform charge distribution with the radius Rq. For 
each £, j, the depth of the central potential has been 
adjusted to the proton separation energy in 17 F, which 
corresponds to a given value of 77. The value of Np(l,j) 
can be extracted from the calculated wave function at 
77 > 100 and it is 10.786, 2.914, and 0.42 for the ls 1/2 , 
Opi/2, and 0d 5 / 2 states, respectively. Figure Q] shows the 

77— dependence of the ratio fiij/Pij. The asymptotic be- 
havior is reached at 77 rj 1 for £ = and 77 rj 10 for 
£ = 2 and the ratio smoothly decreases with 77 suggest- 
ing a polynomial dependence on the proton separation 
energy (or a Maclaurin series in I/77). Guided by this 
result, we write: 



Pnij (rj) = V 1/2 exp(?7 In 77 - 77) ^ 



rj 1 



(16) 



It is worth noting that the term 77 1 ' 2 exp(?7 In 77 — 77), gov- 
erning the rapid growth of SPANC around the threshold, 
is universal for all charged particles (proton, deuteron, a, 
. . .), independent of their quantum state. The depen- 
dence on the structure is contained in the coefficients 
a, of the Maclaurin series, which carry information on 
quantum numbers (n,£,j) of the bound state. In the 
considered example of single-proton states in 17 F, the 




?7 



FIG. 1. The 77— dependence of Ptj/Pij for ls^j OP1/21 an d 
0d 5/ / 2 s.p. proton wave functions in 17 F. 



TABLE I. Coefficients of the fourth-order Maclaurin series 
(116[) fitted to the ratio j3ij//3tj of single-proton WS states in 
F in the range of S p < 10 MeV. For the lsi/2 state, three 
values of diffuseness d (in fm) were used. 







lSl/2 




OP1/2 


04/2 


d 


0.46 


0.58 


0.7 


0.58 


0.58 


ao 


9.934 


10.786 


11.711 


2.914 


0.4200 


(ii 


0.8255 


0.8889 


0.9515 


3.022 


1.185 


(1-2 


5.483 


6.393 


7.444 


1.416 


1.273 


<J3 


-0.9821 


-1.157 


-1.362 


0.1411 


-0.0408 


"4 


0.4598 


0.5635 


0.6945 


0.1676 


0.2035 



coefficients a, were fitted in the range of proton separa- 
tion energies Sp < 10 MeV. An excellent fit has been 
obtained with the first five terms in the expansion (|16l) . 
The resulting values a, (i = 0, . . . , 4) are listed in TablcQ] 
To study the model dependence, for lsi/2 we considered 
three values of WS diffuseness d. In this case, the coef- 
ficients Oj vary by 20-40% if d changes from the value of 
0.46 fm to 0.7 fm. The variations in m with d are further 
reduced if the r.m.s. radius of the potential is kept con- 
stant while changing d. For instance, the changes in ag 
are ~ 1% in the calculations constrained in such a way. 



energy limit of SPANCs for charged 
particles 



2. Large 



Now we consider the limit f(13|) of smaller 77, which 
corresponds to finite separation energies that are small 
enough so that the 77- variations of Mnt can be neglected. 
In this case, the asymptotic part of the wave function 
shows the usual exponential decay: 



U£j(r) ~ faj exp (-K r/r] - 77m(2 K r/r))) 



(17) 



In this case, in order to keep 7V 0X t finite, j3ij has to in- 
crease when 77 decreases. Consequently, when inspect- 
ing fiij as a function of 77 one can expect a minimum 
when 77 changes from small values toward the threshold 

(77 = +00). 



3. Separation energy dependence of SPANCs for bound 
proton wave functions 

We shall now discuss the behavior of j3^j in the full 
energy range. Figure [5] shows Pij{r)) for the p + 16 




FIG. 2. Energy dependence of /3tj for the bound-state proton 
s.p. wave functions in 17 F with IS1/2 ( a )i OP1/2 an d Ops/2 (b); 
and 0d 5 /2 (c). 

capture reaction at E cm = in channels where the pro- 
ton has a relative angular momentum £ — (lsi/2), I = 1 



(Opi/2 and 0p 3 / 2 ), or £ = 2 (0^5/2)- We used the same 
WS potential as in Sec. Ill AT] It is seen that the behavior 
of SPANC as a function of 77 follows general considera- 
tions of Sec. Ill A 21 Namely, with increasing 77, /3ij first 
decreases until a certain minimum value of 77 = ?7 cr it is 
reached, and then increases again as the separation en- 
ergy decreases towards the S p = threshold. We may 
thus conclude that r/ cr it separates the regimes of strong 
(?7 < VcTit) an d weak (77 > ?7 cr it) binding for a given par- 
tial wave and Z a . In general, ?7 cr it scales approximately 
linearly with Z^ and it strongly depends on the angular 
momentum £ of the proton. On the other hand, the de- 
pendence on the channel angular momentum j is weak. 
Fhc magnitude of fyj decreases with £. The case of I = 
and large 77 (and /^=oj=i/2) shown in Fig.^a) is charac- 
teristic of a proton halo. However, as seen in Fig. [2] and 
discussed in Sec. Ill A 21 the large values of (3gj are also 
expected for very bound states, so a large SPANC is not 
an indicator of a proton halo. 

The extreme regimes of SPANC can be characterized 
by the complex turning point z t given by Eq. (JTTJ) . The 
low-binding regime of /3i t j is reached for nr <C \zt\ and 
7/ — > +00. In this region, characterized by the condition 



(77 + ^fTW+l)) T? » KoRf, 



(18) 



as 



the proton wave function behaves asymptotically 
uej(r) oc r 1 ' 4 exp (—2^/2^07-), i.e., its decay is slower than 
exponential. The strong-binding regime is reached at 
nr > \z t \, i.e., 



[v+ y/if 



1) 77<K i? / . 



(19) 



Here, the proton wave function shows the expected ex- 
ponential decay (fT4)) . 



4- Near-threshold behavior of charged particle radiative 
capture cross sections 

Formally, one can discuss the charged particle radia- 
tive capture cross section in the two limits: (i) ?7cm — > 00 

{E CM -> 0), and (ii) 77 -> 00 (S { c a) -> 0). Whereas the 
CM energy of the system (b + c) can be varied experi- 
mentally, the charged particle separation energy is fixed 
for any state of the a-nucleus and, therefore, the limit- 
ing behavior of the radiative capture cross section when 
77 — > 00 cannot be studied experimentally in a single phys- 
ical system. For a fixed value of 77, the radiative capture 
cross-section a cl (j]cm , v) is exponentially reduced in the 
first limit (?7cm — ► 00) as: 



lim 

»7CM->c 



y c7 



Sbc(VCM,v) 



exp(-27r?7CM) 



E, 



CM 



(20) 



As immediately follows from Eq. ((5J) and discussion in 
Sec. Ill ATI in a sub-threshold regime (77 — > 00), the radia- 
tive capture cross-section <j c1 {j)cm-, v) f° r a fixed value of 



r]CM diverges as: 

lim a cl (x Sb c (r]CM,v) « 77exp(277ln(?7) - 2?y). (21) 



rj—too 



The asymptotic behavior of er C7 given by Eqs. (|20p and 
(f2"T|) is expected for any charged particle radiative capture 
reaction. 



5. Survey of experimental proton SPANCs 

Figures [3] and |4] survey the values of fitj and r\ for 
the proton emission channels with t = 1 (p-shell) and 
i = 2 (sd-shell), respectively. The Z a — dependence of 
7/crit is fairly weak in both cases. With the exception of 
the known proton halo 8 B, all other p-shell nuclei belong 
to the class 77 < ?7 cr it- The ground state of 12 N is found 
to have 77 = ?y cr it. The neutron-rich nuclei have very 
small ground-state 77-values; hence, their normalization 
constants j3e = ij(rf) are large (see Fig. (2). An odd-even 
staggering of 77 - due to pairing - leads to an odd-even 
effect in (3i =1 j(Z a ). The staggering in /3tj is stronger for 
nuclei with T z > than in the proton-rich systems with 
T z <0. 
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FIG. 3. Experimental values of 77 (top) and /3ej (bottom) 
for bound proton states with I — 1 and j = 1/2, 3/2 in var- 
ious p-shell nuclei (Z a ,T z ). The curves of 77 cr it and /3 C rit are 
calculated for T z = —1/2 nuclei. See text for details. 

The data for sd-shell nuclei shown in Fig. |4] exhibit 
similar behavior. Most of the particle-stable sd-shell nu- 
clei belong to the class 77 < ?7 cr it- The only system that 



FIG. 4. Similar as in Fig. [3] except for 1 = 2 bound proton 
wave functions in sd-shell nuclei. For certain nuclei, the ra- 
diative proton capture with the channel angular momentum 
j = 3/2, 5/2 cannot populate the ground state of a nucleus a. 
In such cases (encircled symbols), we plot 77 and j3tj for the 
lowest-energy excited state that is populated by the capture 
with £ — 2 protons. 



falls decisively into the regime of weak binding is the 5/2" 
ground state of 23 Al, which has a small separation energ; 
of S„= 1 keV and a very large value of 0e.i (see Ref. [24 



j p — 141 keV and a very large value of Pej 
for a recent discussion). Other nuclei with 77 > ?y cr it are 
29 P, 32 C1, and 36 > 37 K whereas 17 F, 30 S, 33 C1, and 37 Ca 
are situated at the critical line ri CT it(Z a ). A pronounced 
odd-even staggering of 77 and /3gj is seen for both proton- 
rich and neutron-rich systems. The values of SPANCs for 
proton- rich nuclei remain close to /3 C rit- 

The transition from 77 < ?7 cr it to 77 > ?y cr i t can be fur- 
ther explored by studying different particle-stable excited 
states in one nucleus, having different values of 77, popu- 
lated in the capture reaction. Figure EJa) shows exper- 
imental values of 77 as a function of the proton separa- 
tion energy for several excited states of 22 Mg. A transi- 
tion from 77 < 77 cr it to 77 > 77 cr it takes place between the 
J* = 4f and l\ levels. The corresponding SPANCs are 
displayed in Fig.[5][b). Those with £ = 2 show a minimum 



at the l\ state and the same is seen in the £ = channel. 
(One may notice the absence of £ = data points for the 
0^~ and 4+ levels as these states cannot be populated by 
the £ = proton capture.) 
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FIG. 5. Experimental values of r] (a) and /^j (b) for the low- 
energy states of Mg as a function of the proton separation 
energy S p . The values of ?7crit in channels with £ = and 
£ = 2 are indicated. 



B. General properties of SPANCs for neutrons 

The s.p. neutron wave function on the asymptotic re- 
gion is 



Since the external norm A4xt © of ([25]) is finite, the 
near-threshold divergence of J R cxp(—2ftr) dr oc k~ 1 
must be compensated by the K-dependence of SPANC: 



Pe=0j=i/2 ot v 7 ^- 



(26) 



In the limit of the large binding energy, the neutron 
wave function exhibits the exponential behavior: 



utj{r) ~ fitj exp(-w). 



(27) 



Using the reasoning of Sec. Ill A 21 we conclude that fyj 
has to increase with n. Consequently, in the neutron 
case, fiij should monotonically decrease with k all the 
way down to zero. 

To describe fyj for neutrons in a wider energy range, 
one can employ the power series expansion in k: 

oc 

/W/2=5>fco,i/2 K(i+1/2) for^ = 0, (28a) 



i=0 



Pt,3=Yl a tJ K 



«„(i-M) 



for £ ^ 0. (28b) 



i=0 



Figure [5] shows fyj for the bound lsi/2, ®Pi/2, an d 
0c?5/2 neutron s.p. wave functions. The short-dashed lines 
are fits with two leading terms to expansions (|28|). In 

general, the coefficients of- in Eq. (|28|l depend on the 
choice of an average potential. The case £ = is spe- 
cial. Since A4xt -> 1 in the limit S n -> 0, SPANC be- 
comes potential-independent. In particular, the coeffi- 
cient a e 2o 1/2 in Eq- (|2"5)) becomes V2 in the limit of zero 
binding. 



(0) 



For £ = l,j = 1/2, the values of aft are 1.549, 1.612, 
and 1.681 for d equal to 0.46 fm, 0.58 fm, 0.7 fm, respec- 
tively. Variations of a ( ■ with d become insignificant if 
the r.m.s. radius of the potential is kept constant. 



uij(r) ^ Pej (~i £ ) Kr hf(mr), 



(22) 



where h~\ is the spherical Hankel function [22, |23[ . In the 
limit of k — > 

(2/— IV 

H') ^ h+(iK.r)~ { , J" for 1^0. (23) 

[nry 



Using the same arguments as in Sec. Ill A 11 we conclude 
that close to the neutron threshold the neutron SPANC 
behaves as 



/3 £j oc K l for £ ^ 0. 



(24) 



To discuss the special case of £ = 0, we follow the analysis 
of Rcf. j25(. The asymptotic part of the wave function 
for the s-wavc is 



utj(r) ~ pij exp(-Kr). 



(25) 



III. CONTINUUM SHELL MODEL 
DESCRIPTION OF ANCS IN MIRROR SYSTEMS 

An extensive analysis of proton and neutron ANCs for 
light mirror nuclei has been performed using SM and clus- 
ter model wave function and various effective N N inter- 
actions [l|, fl5Hl8j . It has been found that the ratio of 
proton and neutron ANCs for mirror nuclei, 



/v — IC/p/Cjj 



(29) 



is rather insensitive to model details and can be well ap- 
proximated by the expression 



K^K n = 



e **(-«*)F/,_i„(tKpii/) 



K p Rfje(iK n Rf) 



(30) 



where ot(—ir\) is the Coulomb phase shift associated with 
the imaginary momentum ik p , and i^,— «j and jt are the 
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FIG. 6. Ptj for the bound neutron states lsi/2 (a), Opi/2 
(b), and 0d 5 / 2 (c) in 17 as a function of S„. The fits using 
expansions (|28[) are shown by dotted lines. The same WS 
potential as in Fig. [l]was used. 



regular Coulomb and regular Bcssel functions, respec- 
tively. It is to be noted that our expression for IZo differs 
from that of Refs. [l|, [Ty] because of different definitions 
of complex Coulomb wave functions. In our work we fol- 
low the convention of Refs. [2fJ, [27], which results in the 
appearance of an additional factor Q la t(- ir i) m Eq. (|3U|) . 
We shall verify these findings in the CSM. As ANCs 
are impacted by the configuration mixing and contin- 

1/2 

uum coupling through the spectroscopic amplitudes SA 
Q, the OQS framework of the CSM is particularly well- 
suited for the description of the ANCs in mirror nuclei. 
Indeed, since the one-nucleon separation energies in mir- 
ror systems can be appreciably different, the particle con- 
tinuum, both of a resonant and non-resonant character, 
can impact properties of states involved, especially when 
dealing with near-threshold energies. 

The principal difficulty in the formulation of the CSM 
is the treatment of the unbound space of states, i.e., res- 
onances and the non-resonant continuum [lfj [2fj, |28| . 
Therefore, whatever strategy is adopted to formulate 
the configuration interaction approach for OQSs, the key 
points are: (i) the treatment of the s.p. continuum, and 
(ii) the definition of the many-body Fock space. 

Historically, the first approach to formulate the CSM 



in Hilbert space was based on the projection technique 
|29J . Here, the Hilbert space is divided into orthogo- 
nal subspaces Qi cont containing different numbers i con t 
(*cont = 0, 1, • • • ) of particles in the scattering continuum. 
An OQS description of Qo space includes couplings to 
the environment of decay channels through the energy- 
dependent effective Hamiltonian |30 - [33 |: 



W-QoQoiE) = H QoQo + w QoQo( E )i 



(31) 



where Hq q q is the standard SM Hamiltonian describ- 
ing the internal dynamics in the closed quantum system 
approximation and Wq q (E) is the energy-dependent 
continuum coupling term. The Hamiltonian (|31l) is Hcr- 
mitian below the first particle emission threshold and 
complex-symmetric above it. 

The s.p. resonances have to be regularized before in- 
cluding them in the subspace of discrete states. An ef- 
fective way of doing this is based on the anamneses of 
resonances, i.e., images of resonances in the space of 
C 2 — functions |33[. Bound s.p. states, resonance anam- 
neses, and real-energy scattering continuum states form 
a complete s.p. basis J33|- The many-body Hilbert space 
can be represented by Slater determinants spanned by 
this basis. In practical applications of the real-energy 
CSMs, such as SMEC, one uses phcnomcnological argu- 
ments to restrict the number of particles in the scatter- 
ing continuum. Technical difficulties associated with the 
correct treatment of the multiparticle channel wave func- 
tions rapidly grow with i con t; hence, in practical applica- 
tions, the number of particles in the scattering continuum 
has so far not exceeded two [31j ). 

Recently, a different strategy based on the rig ged 
Hilbert space formulation of quantum mechanics [34 1 has 
resulted in the complex-energy CSM (the GSM) |35H37 1 , 
which is a natural generalization of the standard SM for 
unbound systems. In this formulation, the maximum 
number of particles in the scattering continuum is not a 
priori prescribed, but follows from the Schrodinger vari- 
ational principle for the many-body Hamiltonian. The 
s.p. GSM basis is given by the Bcrggrcn ensemble [38j |. 
which consists of Gamow (resonant) states and the non- 
resonant continuum. (For a detailed description of the 
GSM, see Ref. [20[.) The GSM Hamiltonian is Hermi- 
tian. However, since the s.p. vectors have either outgoing 
or scattering asymptotics, the GSM Hamiltonian matrix 
is complex symmetric and its eigenvalues are complex 
above the particle emission threshold. Hence, both real- 
energy and complex-energy CSM formulations lead to a 
non-Hermitian eigenvalue problem above the threshold. 



A. Theoretical framework 

In both GSM and SMEC, we assume that the nucleus 
can be described as a system of A va i = Z va i+N va i valence 
nucleons moving around a closed core A c = Z c + iV c . 



1. GSM framework 

The translationally invariant GSM Hamiltonian, writ- 
ten in intrinsic nucleon-core coordinates of the cluster- 
orbital shell model K59f, can be written as: 



* = £ 



i=\ 






Ui 



E 



V<i + Wp PiPi 



(32) 



where /ii is the reduced mass of either the proton or neu- 
tron (l//i, = 1/m, + l/M c ), Ui is the s.p. potential de- 
scribing the field of the core, V%j is the two-body residual 
interaction between valence nucleons, and the last term 
represents the two-body energy recoil. The particle-core 
interaction is a sum of nuclear and Coulomb terms: 



Ui = U„ 



N 



u: 



c 



(33) 



The nuclear potential Uf is approximated by a Woods- 
Saxon (WS) field with a spin-orbit term [35[, and the 
Coulomb field Up is generated by a Gaussian density of 
Z c core protons [141 ]. Similarly, the residual interaction 
can split into nuclear and Coulomb parts: 



V - - V- N 



Vh 



where V N is the modified surface Gaussian interaction 
(MSGI) [lj| and V c is the two-body Coulomb interaction 
that requires special treatment due to its infinite range 
[Tj,|4(|. Namely, V c is rewritten as: 



ug 



[v 



c u c , 



,HO 



(35) 



where U% _ 1 takes care of the long-range asymptotic 
behavior of the Coulomb interaction. As l/S + US , = 
Uz_ 1: the long-range physics of the Coulomb Hamil- 
tonian is treated almost exactly. The second term in 
Eq. (|35p and the two-body recoil term can be expanded 
in a harmonic oscillator (HO) basis [1J,|41(. In our calcu- 
lations, we took nine HO shells with the oscillator length 
6 = 2 fm. 

The radial overlap integrals were calculated using 
Eq. ^, where the sum over B states runs over the com- 
plete Bcrggren ensemble; hence, the result is independent 
of the s.p. basis representation. The ANC is obtained di- 
rectly from Eq. (Q~J). 

The GSM calculations presented in this paper were 
carried out for p-shell systems 6 Li/ 7 Be (= 6 Li+p — > Be), 
6 Li/ 7 Li (= 6 Li+n -► 7 Li), and 7 Bc/ 8 B, 7 Li/ 8 Li (assum- 
ing 4 Hc core), and sd-shell nuclei 16 0/ 17 F and 16 0/ 17 
(assuming 12 C core). 

As the protons (neutrons) are well bound in 7 Be and 
8 Li ( 7 Li and 8 B), only proton (neutron) bound s.p. shells 
are included in the model space, namely the Op^/2 and 
Opi/2 states. On the other hand, the full p 3 / 2 and 
Pi/2 Bcrggren basis, consisting of resonant and scatter- 
ing states, is taken into account for the neutron (proton) 



space. In this way, completeness is assured in the neu- 
tron (proton) p-space. This is especially important for 
the ground state of the proton halo nucleus 8 B. 

Since the precise description of reaction thresholds is 
crucial for precise determination of radial overlap in- 
tegrals, the parameters defining the WS potential and 
MSGI interaction have been fitted independently for each 
pair of nuclei considered to reproduce the separation en- 
ergy and lowest excited states of the heavier (parent) 
nucleus. The WS potential diffuseness has been assumed 
to be d = 0.65 fm in all cases; WS radius has been set 
to Ro = 2fm for A = 6 — 8 nuclei and i?o = 2.8 fm for 
A = 16/17 pairs; and spin-orbit strength was taken as 7.5 
McV for A = 6 - 8 nuclei, and 7.920 MeV (proton) and 
8.463 MeV (neutron) for A = 16/17 pairs. For A = 6/7 
pairs, the WS strength is 45.455 MeV for 6 Li/ 7 Be (pro- 
tons) and 6 Li/ 7 Li (neutrons), and 55 MeV for 6 Li/ 7 Li 
(protons) and 6 Li/ 7 Be (neutrons). For A = 7/8 pairs, the 
corresponding values are 46.273 MeV (protons) and 65 
MeV (neutrons) for 7 Be/ 8 B, and 75 MeV (protons) and 
45.799 MeV (neutrons) for 7 Li/ 8 Li. For A = 16/17 pairs, 
the WS strength is 46.427 MeV and 46.034 MeV for pro- 
tons and neutrons, respectively. The coupling strengths 
of MSGI depend on the quantum numbers J and T of 
the nucleon pair 



141. The T = constants V 



N 
J,T=0 



are 



(34) listed in Table |TT] and the T = 1 constants are given by 



V N 
V J,T= 



0.5Vj% 



TABLE II. Coupling constants Vf T=0 of the MSGI residual 
interaction (in MeV fm 3 ) [1J] for the considered pairs of nuclei 
and all allowed values of J. The J— 4 and 5 coupling constants 
in A — 16/17 systems have been set to —9.0. 





J = 


J=l 


J = 2 


J = 3 


6 Li/ 7 Be 


-34.000 


-23.200 


-24.300 


-17.384 


6 Li/ 7 Li 


-34.000 


-23.200 


-24.000 


-17.000 


7 Be/ 8 B 


-25.593 


-16.007 


-10.288 


-12.400 


7 Li/ 8 Li 


-24.662 


-14.957 


-10.592 


-12.400 


16 0/ 17 F 


-16.540 


-12.564 


-5.870 


-9.000 


16Q/17Q 


-16.540 


-12.564 


-5.870 


-9.000 



For A = 16/17 systems, the model space is the same for 
protons and neutrons: it consists of the resonant shells 
0pi/2, 0^/2 > and lsi/2, and the S]y 2 and d 5 / 2 scattering 
continua. The p\/2 and d%/2 partial waves are discarded 
as they do not impact the asymptotic behavior of the 0J 
ground state of 16 and the 5/2^" and 1/2] 1 " states of 17 
and 17 F. In the pole approximation, in which only reso- 
nant states are considered, amplitudes of 2p-2h and 4p-4h 



Op 



1/2 



sd excitations are of the order of 10 and 10 



respectively. Consequently, as configuration mixing ef- 
fects in the states considered arc weak, only two particle- 
two hole excitations have been allowed from the proton 
and neutron Opi/2 states, and only one particle has been 
allowed to occupy S1/2 and d 5 i 2 scattering states. In all 
cases, scattering contours have been discrctized utilizing 
a Gauss-Lcgendrc quadrature. We have checked that a 



45-point discretization for A = 6/7 (and 30-point dis- 
cretization for other systems) is sufficient to maintain a 
correct asymptotic behavior oi Ig c . e j(r) up to at least ~ 8 
fm. For the s.p. basis, we took a Gamow Hartree-Fock 
(GHF) ensemble [37| corresponding to a parent nucleus. 
Sphericity of the GHF potential is guaranteed by the use 
of the uniform filling approximation. By taking the GHF 
basis we minimize configuration mixing. The many-body 
GSM states have been determined by a diagonalization of 
the GSM matrix using the Davidson method extended to 
complex-valued Hamiltonians. The identification of the 
outgoing GSM states has been carried out by applying 
the overlap method |36| ■ 



TABLE III. Calculated and experimental separation energies 
for the considered pairs of nuclei. 





r{A- 


1) 


r(A) 


£gsm (MeV) 


£ oxp (MeV) 


6 Li/ 7 Be 


2+ 




3/2~ 


5.606 


5.606 


6 Li/ 7 Li 


2+ 




3/2- 


7.250 


7.250 


7 Be/ 8 B 


3/2" 




2+ 


0.137 


0.137 


7 Li/ 8 Li 


3/2~ 




2 + 


2.034 


2.033 


16 0/ 17 F 


0+ 




5/2+ 


0.601 


0.600 


1G 0/ 17q 


0+ 




5/2+ 


4.143 


4.143 



The parameter optimization was carried out using the 
multidimensional Newton method. The fine-tuning was 
done manually to adjust the thresholds to experimental 
value with a 1 keV precision. For the excited states, we 
reproduce experimental data with a precision of a few 
tens of keV. The calculated separation energies are com- 
pared to experiment in Table Hill while Table HVl displays 
excitation energies and widths of the lowest excited states 
in the nuclei considered. 



TABLE IV. Calculated and experimental excitation energies 
(in MeV) and widths (in keV) of the first excited states of the 
parent nucleus. The approximate width Iqq^ of Eq. (|47[) is 
shown for comparison. 





r 


-Egsm 


-C'cxp 


Tgsm 


p(app) 

1 GSM 


r 

i cxp 


7 Be 


1/2- 


0.439 


0.429 











7 Li 


1/2- 


0.459 


0.478 











8 B 


1+ 


0.771 


0.770 


24 


23 


36 


8 B 


3+ 


2.278 


2.320 


275 


258 


350 


8 Li 


1+ 


0.992 


0.981 











8 Li 


3+ 


2.222 


2.255 


21 


21 


32 


iv F 


1/2+ 


0.495 


0.495 











u O 


1/2+ 


0.870 


0.871 












2. SMEC framework 

In this study, the scattering environment is provided by 
one-nucleon decay channels, i.e. we solve the Schrddingcr 



equation in the function space Qq © Q\. The SMEC 
solutions in Qq are found by solving the eigenproblem 
for the non-Hcrmitian Hamiltonian (J31I) : 



UQ oQo \y a ) = £ a {E,V )\* a ) 

(* a \H QoQo =£* a (E 7 V )(*a\ 



(36) 



in the biorthogonal basis: (^fa\^p) = Safi- As usual, 
left I^q) and right |\l/ s ) eigenvectors are related by the 
complex conjugation. In Eq. (|36[) . E and Vq stand for a 
scattering energy and a (real) continuum coupling con- 
stant in the coupling terms 'H.Q ay Q 1 and 'Hq 1 q . 



W Q Qo( E ) - H QoQi G o, ( E ) H QiQo 



M, 



(37) 



-.(+), 



where Gq (E) is the one-nucleon Green's function. The 
energy scale is defined by the position of the one-nucleon 
emission threshold. At resonance, the eigenvalue of the 
effective Hamiltonian can be identified with the narrow 
pole of the scattering matrix (the S'-matrix). 

Inside of the interaction region, dominant contribu- 
tions to the full solution of the Schrddinger equation in 
Qo © Qi are given by the eigenfunctions of T-Lq q (E). 
This is the main reason why eigenfunctions of the non- 
Hcrmitian effective Hamiltonian are essential to under- 
stand properties of the OQS. The SMEC eigenvectors 
^f a are related to the eigenstates $j of the closed quan- 
tum system Hamiltonian Hq q by a linear orthogonal 
transformation: 



*a = 5> 



■CLJ^j- 



(38) 



In our SMEC calculations, for the effective SM Hamil- 
tonian Hq q we took the Cohen-Kurath interaction [42| 
for A = 6 — 8 systems and the ZBM effective interaction 
43] for sef-shcll nuclei. The continuum-coupling term 
(|3"T)) was approximated by means of the Wigner-Bartlett 
contact interaction: 



V 12 =V [a+{l- a)P^ 2 ] 8 (n - r 2 ) 



(39) 



where P" 2 is the spin exchange operator and a = 0.73. 
The magnitude of the continuum coupling varies depend- 
ing on the structure of SM wave function in a target 
nucleus with (A — 1) nucleons and the energy of the low- 
est one-nucleon emission threshold, which is fixed at the 
experimental value in all calculations. SMEC is partic- 
ularly suited for studies of the qualitative effects of the 
continuum coupling because the strength of this coupling 
can be changed continuously from a SM limit (Vq = 0) 
to physically relevant values. 

The expectation value of any operator O can be calcu- 
lated as: 

(6) = (tf 6 |6|* Q ) , (40) 

In case of the spectroscopic factor one has: 

6 = a)\t)(t\ ai (41) 
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where |t) is the target state of the (A — l)-system. For a 
single SM configuration, the ANC is proportional to the 
square root of the spectroscopic factor ([7|). In SMEC, the 
spectroscopic factors depend on the total energy E of the 
system and exhibit characteristic near-threshold varia- 
tions that depend on the transferred angular momentum. 
In the multichannel representation of a many-body sys- 
tem, the flux conservation imposes an intricate interde- 
pendence between various spectroscopic factors not only 
on E but also on the strength of the continuum coupling 
Vq . This salient dependence of the ANC on the contin- 
uum coupling strength is a quantal effect, beyond the 
generic features discussed in Sec. HU 

In the case of multi-channel coupling, the squared 
norm ([3]) becomes 



iq 2 = £k%| 2 &. 



(42) 



t-j 



The ratio of proton and neutron ANCs for mirror nu- 
clei (P5|) can be directly computed by means of Eq. (|4"2")l 
applied to charged-particle and neutral-particle radiative 
capture reactions. 



B. CSM description of ANCs for bound states 

In this section, we discuss the ANCs corresponding to 
single-nucleon capture reactions between bound states. 
We first present our GSM results (E C m = 0). The ANCs 
in GSM can be directly extracted from the calculated 
radial overlap integrals by fitting their tail to Whittaker 
functions at large values of r (=7-8 fm). We can use such 
a direct method of extraction because the asymptotic 
behavior of Ibc-ej( r ) * s weu controlled in GSM. 

Figure [7] shows how this procedure works for the ra- 
dial overlap integrals corresponding to t = 1 (0pi/2 an d 
OP3/2) protons and neutrons in the 2 + ground states of 
the mirror nuclei 8 B, 8 Li. This example is non trivial as 
the configuration mixing is appreciable, and the g.s. of 
8 _B is a proton halo - as seen from the extended tail of 
the overlap function in Fig.fTJa). 

The second example presented in Fig. [8] pertains to 
the first excited J 71 " = l/2f subthreshold halo state in 
17 F and 17 0, as seen in Table IIVI Here, both proton 
and neutron overlap functions are very extended. Never- 
theless, the extraction of the ANCs does not cause any 
problems. 

The GSM predictions for ANC are listed in Table E] 
together with the Variational Monte Carlo (VMC) re- 
sults of Ref. Q and experimental data. One can see that 
the values of ANC for the P1/2 proton and neutron par- 
tial waves are very different in GSM and VMC. While 
VMC values of ANCs in p-shell nuclei are usually closer 
to experiment for individual partial waves, both GSM 
and VMC perform well when the norm ([3]) is consid- 
ered. The GSM results for the 5/2 + partial wave in the 
A = 16 — > 17 capture are very close to the data. This 
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r (fm) 

FIG. 7. Radial overlap integrals I^ c; ij(r) calculated in GSM 
for 7 Be 3 / 2 - +p — > 8 B 2 + (a) and 7 Li 3 / 2 - + n — > 8 Li 2 + (b) for 
I — 1 and j = 3/2. The tail of the radial overlap integral 
is fitted by the Whittaker function (dashed line) to extract 

ANC. 
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and 16 O 0+ + n -)• 17 Oi /2 + (b) for I = 0. 



does not come as a surprise as the associated spectro- 
scopic factors are almost equal to one, i.e., the ANC and 
SPANC values are practically identical. 

The ratios K (EH) predicted in GSM and VMC are 
displayed in Table IVII They are compared to experimen- 
tal data and the approximate expression TZq of Eq. (|30[) . 



11 



TABLE V. GSM predictions for ANCs in mirror systems com- 
pared to VMC results [f| and experimental data p], l44T - |47l | . 



For b Li 1+ + n 



'Li 



3/2" 



two experimental values have been 



reported [4J] . For A — 6 — 8 systems, we show individual con- 
tributions from P1/2 and P3/2 channels and their Hermitian 
norm ((3]). 



Overlap 


f 


GSM 


VMC 


Exp. 




1/2- 


0.431 


1.870 


— 


6 Li 1+ -* 7 Be 3/2 _ 


3/2 _ 


1.499 


2.150 


— 




1/2- +3/2- 


1.559 


2.850 


— 




1/2- 


0.422 


1.652 


— 


6t ■ "71: 
L>1 1 + — > L1I3/2- 


3/2- 

1/2-+3/2 - 


1.456 
1.559 


1.890 
2.510 


1.86 ± 0.06 
2.57 ± 0.06 




1/2- 


0.049 


0.246 


0.23 ± 0.01 


Be 3/2 - — > B 2 + 


3/2- 


0.765 


0.691 


0.64 ± 0.03 




1/2-+3/2 - 


0.767 


0.733 


0.68 ± 0.04 




1/2- 


0.041 


0.218 


0.22 ± 0.01 


7 T ■ n 8t: 
1j1 3/2- — ' ijl 2+ 


3/2- 


0.750 


0.618 


0.62 ± 0.03 




1/2- +3/2- 


0.752 


0.655 


0.66 ± 0.03 


16 n P. 17-p 


5/2+ 


0.880 


— 


0.95 ± 0.09 


1/2+ 


73.74 


— 


— 




5/2+ 


0.805 


— 


0.82 ± 0.01 


D o 0+ -> G> 


1/2+ 


2.785 


— 


— 



TABLE VI. Ratio H ([29]) calculated in GSM and VMC 0, 
compared with estimate 7?-o (|30|) and experiment for the mir- 
ror systems of Table [V] 



Mirror pair 


f 


75.GSM 


Kvmc 


Ko 


7?-Exp. 




1/2- 


1.04 


1.28 


1.06 


— 


( 7 Be/ 7 Li) 3/2 - 


3/2- 


1.06 


1.29 


1.06 


— 




1/2-+3/2 - 


1.06 


1.29 


1.06 


— 




1/2- 


1.39 


1.27 


1.12 


1.08 ± 0.18 


( 8 Be/ 8 Li) 2+ 


3/2- 


1.04 


1.25 


1.12 


1.08 ± 0.15 




1/2-+3/2" 


1.04 


1.25 


1.12 


1.08 ± 0.15 


( 17 F/ 17 0)^ 


5/2+ 
1/2+ 


1.20 
701 


— 


1.22 
796 


1.33 ± 0.20 



(The predictions of microscopic cluster models can be 
found in Refs. [l[ and jla|.) As discussed in Rcf. [l|, as 
compared to ANCs, values of 1Z exhibit less model de- 
pendence. Overall, predicted ratios 71 are fairly close to 
the estimate TZq and experiment. 

In the following, we shall discuss the mirror ANCs us- 
ing SMEC. The advantage of the projection technique 
used in SMEC is that the continuum coupling can be 
switched off and, therefore, the effect of the environment 
of decay channels on the mirror ANCs can be studied 
separately from the effect of internal dynamics governed 
by the effective interaction. The SMEC calculations pre- 
sented in this paper were carried out for p- and sd- 



shell bound-state mirror reactions: 6 Li/ 7 Be and 6 Li/ 7 Li; 
7 Be/ 8 B and 7 Li/ 8 Li; n C/ 12 N and U B/ 12 B; 16 0/ 17 F and 
16 0/ 17 0; and 17 F/ 18 Ne and 17 0/ 18 0. 



TABLE VII. SMEC predictions for ANCs in mirror systems 
calculated with Cohen-Kurath (A = 7, 8, 12) and ZBM (A = 
17, 18) interactions for three values of the continuum coupling 
strength Vo (in GeVfm 3 ): (SM limit), -0.65, and -1.30, 
all in the physically relevant range of continuum-coupling. 



Ov 


erlap 


r 


(SM) 


0.65 


1.30 


6 Li 1+ ^ 


» 7 Be 3/2 - 


1/2- 

3/2" 
1/2- +3/2- 


1.017 
1.279 
1.663 


1.024 
1.282 
1.641 


1.046 
1.291 
1.661 


6 Li 1+ - 


"+ ljl 3/2- 


1/2- 
3/2" 

l/2~+3/2^ 


0.988 
1.243 

1.588 


0.995 
1.246 
1.595 


1.016 
1.254 
1.614 


7 Be 3/2 - 


— > h> 2 + 


1/2- 

3/2- 

1/2- +3/2- 


0.164 
0.702 
0.721 


0.149 
0.702 

0.717 


0.088 
0.700 
0.705 


L13/2- 


— > Ll2+ 


1/2- 

3/2" 
l/2~+3/2" 


0.160 
0.685 
0.704 


0.151 
0.685 
0.702 


0.115 
0.684 
0.704 


^3/2- 


— > iN 1 + 


1/2- 

3/2" 

1/2- +3/2 _ 


1.122 
0.535 
1.242 


1.126 
0.530 
1.244 


1.127 
0.527 
1.244 


lln 

^3/2- 


n , 12 R 
— ¥ B 1+ 


1/2- 

3/2- 

1/2" +3/2" 


0.971 
0.465 
1.077 


0.975 
0.459 
1.078 


0.977 
0.454 
1.077 


16 n 


A 17 F-x 


5/2+ 
1/2+ 


0.869 
75.91 


0.867 
76.01 


0.863 
76.09 


16 n 


A 17 0,. 


5/2+ 
1/2+ 


0.788 
2.771 


0.787 
2.773 


0.783 
2.776 


17 p 

*5/2+ 


A 18 Ne 2+ 


1/2+ 

5/2+ 

1/2++5/2+ 


4.500 
1.739 
4.823 


4.610 
1.758 
4.934 


4.961 
1.702 
5.245 


17 o 5/2+ 


A 18 o 2+ 


1/2+ 

5/2+ 

1/2++5/2+ 


2.801 
1.504 
3.174 


2.796 
1.544 
3.194 


2.794 
1.580 
3.210 


^5/2+ 


A 18 Ne 2+ 


1/2+ 

5/2+ 

1/2++5/2+ 


13.36 
1.644 
13.46 


13.92 

1.820 
14.04 


12.36 
2.299 

12.57 


O5/2+ 


A 18 o 2+ 


1/2+ 

5/2+ 

1/2++5/2+ 


2.492 
0.705 
2.590 


2.674 
0.683 
2.760 


2.774 
0.679 
2.856 



Table IVIII contains the summary of SMEC predictions 
for ANCs in p- and sd-shell mirror nuclei. To illustrate 
the impact of continuum coupling, we varied the contin- 
uum coupling strength Vo in the physically relevant range 
from (the SM limit) to - 1.30 G eV fm 3 . The resulting 
ratios TZ are listed in Table I Villi 

It is seen that the dependence of ANCs and TZ on 
the continuum coupling is usually very weak. Indeed, 
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TABLE VIII. Ratio H 
systems of Table NTH 



calculated in SMEC for the mirror 



Mirror pair 


r 


7?-SM 


7^-0.65 


7^1.30 




1/2- 


1.059 


1.059 


1.060 


( 7 Be/ 7 Li) 3/2 _ 


3/2 - 


1.058 


1.058 


1.059 




1/2- +3/2- 


1.058 


1.058 


1.060 




1/2- 


1.055 


0.974 


0.584 


( 8 Be/ 8 Li) 2+ 


3/2- 


1.048 


1.048 


1.047 




1/2 - +3/2- 


1.049 


1.045 


1.034 




1/2- 


1.333 


1.333 


1.331 


( 12 N/ 12 B) 1+ 


3/2- 


1.323 


1.333 


1.347 




1/2- +3/2- 


1.331 


1.333 


1.334 


( 17 F/ 17 0)^ 


5/2+ 
1/2+ 


1.216 
750.7 


1.215 
751.4 


1.213 
751.3 




1/2+ 


2.580 


2.719 


3.153 


( 18 Ne/ ls O) 2 + 


5/2+ 


1.336 


1.295 


1.160 




1/2++5/2+ 


2.302 


2.386 


2.670 




1/2+ 


28.74 


27.10 


19.85 


( 18 Ne/ ls O) 22+ 


5/2+ 


5.438 


7.099 


11.46 




1/2++5/2+ 


27.02 


25.88 


19.37 



in most considered cases the effect of the continuum cou- 
pling on ANCs does not exceed a few percent, and it is 
even smaller -a few per mil - for 1Z. It is instructive to 
compare ANCs of Tables [V] and IVII| and the ratios 1Z of 
Tables IVII and IVIIII for the mirror pairs ( 7 Be/ 7 Li) 3 / 2 - , 
( 8 Be/ 8 Li) 2 +, and ( 17 F/ 17 0)^. The GSM and SMEC re- 
sults are extremely consistent when it comes to the the 
total ANCs (|42]) and their ratios. 

As will be shown below, the effect of the continuum 
mixing depends on the distribution of spectroscopic fac- 
tors in SM states coupled to the decay channel. This 
distribution is shown in Fig. [9] for selected examples dis- 
cussed in this section. 

A typical example, illustrated in Fig. [TTJ1 shows 
the variation of ANCs with Vb for the mirror pair 
( 7 Be/ 7 Li) 3 /2- and the mirror capture reactions: 6 Li!+ + 



P 



'Be 



3/2" 



7 and Lii+ 



Li 3 /2- + 7- Tne 



change of C n e. with Vb is due to the mixing of different 
3/2~ SM states caused by the continuum coupling. This 
external mixing of SM states changes the spectroscopic 
amplitudes S Pl/2 and S P3/2 in C Pl/2 and C P3/2 ANCs, re- 
spectively. As seen in Fig.[9ja), the spectroscopic factors 
Sq Pi/2 and <Sop 3/2 m J* — 3/2~ states (i = 1, 2, 3) have all 
comparable values. The continuum mixing in this case is 
small, on the order of 2% (see Table |VlTJ). 

Table IVIIII shows that in spite of different proton 
and neutron separation energies in the mirror pair 
( 7 Be/ 7 Li) 3 / 2 - , the effect of the continuum coupling on 
the ratio 1Z for the 1/2", 3/2" partial waves, and the 
squared norm (|4"2")l is exceedingly small. 

We now consider the mirror pair ( 8 B/ 8 Li) 2 + - in- 
volving a proton halo 8 B - and the mirror reactions: 



¥ ° 

.<& 1 



&-2 







( a ) ( 7 Be/ 8 Li) 3/r Or (17 F/ 17 ) (C) 




o 

b 

o 

& 1 
-2 -2 



shell model state 

FIG. 9. (Color online) Distribution of spectroscopic strength 
for SM states coupled to the same decay channel in mirror 
pairs: ( 7 Be/ 7 Li) 3/2 - (a); ( 8 B/ 8 Li) 2+ (b); ( 17 F/ 17 0) 1/2+ (c); 
and ( 18 Ne/ 18 0) 2+ (d). 
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1.26 
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-1000 -500 

V (MeV fm 3 ) 

FIG 10. ANCs for mirror nuclei 7 Be 3 / 2 - (solid line) and 
7 Li 3 / 2 - (dashed line) calculated in SMEC as a function of 
the continuum-coupling strength Vb- Since the target nucleus 
6 Li has J x = 1+, two partial waves are possible: pi/ 2 (a) and 
P3/2 (b). 



'Be 



3/2- 



% 



7 and 7 Li ; 



3/2- 



3 Li 2 + + 7. 



The target nuclei have J T = 3/2~; hence, they can be 
coupled to the final 2+ state through pi/ 2 or p 3 / 2 waves. 
Figure [TT] shows the corresponding ANCs. The effect of 
the continuum coupling is rather important for a small 
component P1/2 and practically negligible in p 3 / 2 . No- 
tice also a rather strong - and different - dependence of 
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FIG. 12. Similar as in Fig. [10] but for mirror nuclei Ne + 

(solid line) and 18 2 + (dashed line). Since the target A = 17 

nucleus has J^ = 5/2 + , two partial waves are possible: d 5 / 2 
(left panels) and S\i% (right panels). 



FIG. 11. Similar as in Fig. [10] but for mirror nuclei 8 B 2 + 
(solid line) and 8 Li 2 + (dashed line). Since the target A — 1 
nucleus has ,P — 3/2~, two partial waves are possible: p\/i 
(a) andp 3/2 (b). 



C Pl/2 on Vq in mirror systems. This different response to 
the continuum-coupling can be traced back to a differ- 
ent distribution of SM spectroscopic factors in the three 
lowest 2 + states; see Fig. |9][b). As the spectroscopic 
factor So P3/2 0I the 2f state is close to 1, the state 2+ 
is aligned with the decay channel already at Vo=0 and 
no further redistribution of spectroscopic strength is pos- 
sible through the continuum coupling. The situation is 



different for S 



OPl/2 



In this case, the second 2 + state has 



the largest spectroscopic factor and the external mixing 
leads to a redistribution of spectroscopic strength; hence, 
a change in ANC. The ratio of ANCs for the mirror pair 
( 8 B/ 8 Li) 2+ is shown in Table NTU\ The variation of K 
with Vo is of the order of 1 percent, and practically the 
whole effect is due to the "px/i wave. 

As a third example, relevant in the context of GSM 
analysis, we shall consider the mirror pair ( 17 F/ 17 0) in 
J n = 5/2+ ground state and in the first excited proton 
halo state l/2 + . In the J v = 5/2 + ground state, both 
17 F and 17 couple to 16 through the d 5 / 2 wave. The 
0(^5/2 spectroscopic strength is practically localized in the 

lowest 5/2 state, i.e., the ground state is aligned with 
a decay channel. The same is true for the excited 1/2+ 
state, which practically exhausts the lsx/2 spectroscopic 
strength, see Fig. [He). Consequently, as seen in Tables 
IVIII and IVIII1 the continuum coupling is negligible in the 
( 17 F/ 17 0) case. 

In the following, we shall discuss the mirror pair 
( 18 Ne/ 18 0) in the two lowest 2 + states, where the contin- 
uum coupling impacts ANCs significantly. In this case, 



the corresponding mirror reactions are: 



5/2 



+ P 



IN 



NCn 



7 and O 



i7 n + 



is 



5/2 



O2 + 7> where i = 1,2. 
In 18 0, both 2 + states are well bound, whereas in 18 Ne 
the state l\ is close to the proton threshold. It has been 
shown [48[ that this state aligns strongly with the decay 
channel due to the continuum mixing of different 2 + SM 
states. 

Figurc[T2"lshows the mirror ANCs for d 5 / 2 and Si/ 2 par- 
tial waves. It is interesting to notice that with increas- 
ing continuum coupling, Cd 5/2 first increases and then 
strongly decreases in 18 Nc, whereas it steadily increases 
in 18 0. As seen in Table Ivm the overall variations in 
Cd 5/2 are ~ 8%. Even stronger variations with Vq are 
seen for C Sl/2 . In the studied range of Vo values, C Sl/2 
changes by almost ~ 10% in 18 Ne while it varies by ~ 1% 
in 18 0. A different behavior of Cd 5/2 and C Sl/2 results in 
a particularly strong variation (~ 15%) of TZ for the 2f 
state; sec Table Ivmi This behavior can be attributed to 
the distribution of spectroscopic strength 0^5/2 (lsi/2), 
which is primarily concentrated in the two (three) lowest 
2+ SM states. 

Figure Q2] shows the ANCs for d 5 / 2 and S1/2 partial 
waves for the 2\ state in 18 Ne and 18 0. Also in this case, 
ANCs are strongly affected by the continuum coupling: 
the ratio 1Z changes by almost 50% in the considered 
range of Vq , as seen in Table IVIII1 The distribution of 
spectroscopic strength 0d 5 / 2 (lsi/2) in this case is con- 
centrated in the two (three) lowest 2 + SM states; see Fig. 
Hd). 

The realistic examples of SMEC calculations of ANCs 
presented in this section demonstrate that the distribu- 
tion of spectroscopic strength over an ensemble of J T SM 
states is crucial for determining the continuum coupling 
effect on ANCs. If the spectroscopic strength is strongly 
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FIG. 13. Similar as in Fig. ll2l but for J n = 2% states in mirror 
nuclei. 



localized in one state, as in 17 F and 17 0, or very broadly 
distributed, as in 7 Be and 7 Li, then the corresponding 
ANC is fairly insensitive to the continuum coupling. On 
the other hand, if the spectroscopic strength is concen- 
trated in several close-lying SM states, like in 18 Ne and 
18 O, both mirror ANCs and their ratios may strongly 
depend on the coupling to the continuum - in particular 
if the state of interest lies close to the particle-emission 
threshold. 

The distribution of spectroscopic strength strongly de- 
pends on the effective nucleon-nucleon interaction. In 
this sense, the quantitative effect of the continuum cou- 
pling on ANCs is strongly interaction-dependent. One 
should keep this in mind when making predictions about 
mirror reaction cross-sections. For a given model space 
and SM interaction, the relative importance of the con- 
tinuum coupling on mirror ANCs can be a priori assessed 
by calculating spectroscopic amplitudes and their distri- 
bution in a standard SM. However, the effects of contin- 
uum coupling cannot be considered in an isolation from 
the optimization of the SM interaction to spectroscopic 
and reaction observables within a unified framework. By 
doing so, an inherent arbitrariness associated with pre- 
dictions of ANCs can be reduced. 



C. CSM description of ANCs for unbound states 



The definition of ANCs via Eqs. (|6I7[) is no longer ap- 
propriate for negative separation energies, i.e., when the 
state of a nucleus a (A-particle system) is unbound with 
respect to the nucleus b (A — 1-particle system). Indeed, 
in this case k becomes complex and the Whittaker func- 
tion becomes complex as well. The imaginary part of the 
Whittaker function is not vanishing even at the limit of 
vanishing width and the associated ANCs are complex. 



1. ANC of a complex- energy state and its relation to the 
particle width 

A suitable definition of the ANC for positive sep- 
aration energies involves the outgoing Coulomb wave 
function Hf(kr), where k = (—2/iSa/h 2 ) 1 / 2 and rj = 
Zt,Z c e 2 fi/h 2 k: 



1 



Iwfr) ~ -CijH^kr). 



(43) 



Note that k and r\ are complex [T2|, as the state in A- 
particle systems is unbound. At the limit of vanishing 
width, fc, Ifr..ij and Hf(kr) become real, so that CV, 
given by Eq. (|4"3")l becomes real as well. 

For narrow resonances, ANCs can be related to the 
particle width B. However, the derivation of this rela- 
tionship in Ref. [4[ relies on the R-matrix theory - not 
used in the context of GSM - so it is useful to recall 
the derivation for the Gamow states. The overlap func- 
tion I£ c .(j(r) defined in Eq. ([1]) obeys a Schrodinger-like 
equation, albeit inhomogencous [49|. However, separat- 
ing the full interaction into a one-body term and a two- 
body residual interaction, the source term can be decom- 
posed into a dominant homogeneous part and a residual 
inhomogeneous part, and the latter can be absorbed into 
the homogeneous part. This approximation has been 
tested successfully in Ref. 12( for both bound and un- 
bound states. Moreover, as only narrow resonant states 
are involved, we will consider that the potential entering 
the equation defining I§ c .£j(r) is real. This simplifica- 
tion breaks down for resonant states bearing a sizeable 
width, for which complex potentials must be used [12| . 
but is sound for narrow resonances. 

Under these assumptions, one can easily derive the re- 
lation between ANC and partial width [50| ■ The function 
I(r) = rl£ c . £i (r) is a solution of the Schrodinger equation: 



I"{r) 



e{e + i) 



v(r) 



I(r), 



(44) 



where v(r) is real and local. The boundary conditions 
defining I(r) are I(r = 0) = and Eq. (|43| at large r. 
The continuity equation for I(r) implies that 

I*(r)i'(r)-i'*{r)I(r) 

(45) 



(jfc* -AT) / \I{s)\ z ds. 



Taking r in the asymptotic zone in which (|43|) applies, 
noticing that k* — k 2 is proportional to the partial 
width T^j, and utilizing the standard mirror relation 
for Coulomb wave functions Hf{z) = Hj ' ,(z*), (both 
functions obey the same differential equation and behave 
as exp(— iz* + irf ln(2z*)) for \z\ — > +oo), one obtains: 



r <y - 



kH~ r (k*r)H+ ri (kr)' - k*H~ rf (k*r)' ' H+Jkr) 



2if \I{s)\ 2 ds 



h 2 
x—\C e 



(46) 
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where fi is the effective mass of the particle. To get rid of 
the explicit r-dependence in Eq. (|4"u]) . further approxima- 
tions are necessary [51|. Neglecting ^s(k) in the Coulomb 
wave functions of the numerator of Eq. (|46|) implies that 
their Wronskian becomes equal to 2i3t(fc). Moreover, as 
I(r) has a quasi-bound state character, it decreases ex- 
ponentially along the real r-axis (unless r becomes ex- 
tremely large, which we do not consider here), so that 
the integral in the denominator is almost equal to one 
when r is chosen in the asymptotic region. Under these 
assumptions, valid for narrow resonances, Eq. (|46[) sim- 
plifies to: 



|Czi| 2 3l(fc), 



(47) 



which is the same expression as that obtained in Ref . [J] , 
even though approximations and boundary conditions 
are different in the real-energy R-matrix approach and 
complex-energy Gamow-state formalism |5l| . 

Expressing the total width T in terms of the sum of 
partial widths Tij gives total ANC constant C ([3]): 



C 




(48) 



2. CSM description of ANCs for unbound states 

Figure [14] compares radial overlap integrals <j2j) for the 

excited J 71 ' = if states of 8 B and 8 Li in the p^/2 chan- 
nel. Since the if state of 8 B is a narrow one-proton 
resonance, the proton overlap integral acquires a small 
imaginary part. Figure [15] shows the radial overlap in- 
tegrals for a 3+ broad resonance in 8 B and a narrow 
mirror resonance in 8 Li. The tails of real and imaginary 
parts of radial overlap integrals are fitted with the out- 
going Coulomb wave functions of a complex argument k. 
Table IIXI displays proton and neutron ANCs and their 
ratio 1Z for the first excited Jf = if state in 8 B and 8 Li. 



TABLE IX. Proton and neutron ANCs, 11, and TZv (EH calcu- 
lated in GSM for the J 71 = 1^ excited state in S B (one-proton 
resonance) and 8 Li. See text for details. 



J 



1/2" 
3/2 _ 

1/2-+3/2" 



C n 



C n 



7£gsm Ti-r 



0.0322-i0.00138 0.1379 

0.0442-i0.00273 0.2090 

0.0547 0.2504 



0.0545 0.0021 
0.0449 0.0018 
0.0478 0.0019 



3. SMEC description of ANCs for unbound states 

To relate ANCs in bound-unbound mirror pairs and 
extract the proton decay width from the neutron ANC in 
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FIG. 14. Similar as in Fig. [7] but for 7 Be 3/2 - +p -> 8 B 1+ (a) 
and 7 Li 3 , 2 - + n — > 8 Li 1+ (b). 
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FIG 15. Similar as in Fig. [I4]but for the J n = 3^ resonance. 
In 8 B this state is a broad one-proton resonance which results 
in a large imaginary part of radial overlap integral. 

the bound mirror state, Timofeyuk et al. [l| introduced 
a quantity: 



n r 



i r r 



(49) 



he |C„| 2 ' 
where \C n \ is given by Eq. (|42|) and T p is the total width: 

r p = J2 r a\Sij\, (50) 



16 



where r^j is the partial proton width of a Gamow state 
with an quantum numbers I and j. The modulus of a 
spectroscopic factor Sgj is taken in order to ensure that 
T p remains positive after the coupling to the particle con- 
tinuum. For narrow resonances, using Eq. (|47|l. one can 
express IZr in terms of proton and neutron ANCs: 



n r = \c p /c n \ 2 R(*p) 



he 

flpC 2 



(51) 



As an example, let us consider the previously discussed 
case of the proton 1 + resonance in 8 B and its bound mir- 
ror analog in 8 Li. Figure [12] shows the continuum cou- 
pling strength dependence of IZr computed in SMEC. 
The SMEC prediction is compared with the approximate 
formula (|50[) using the spectroscopic factors obtained in 
SMEC. The shaded region correspond to the estimate of 
Ref. [l[ based on experimental data for mirror states in 
8 B and 8 Li. The SM prediction is ftp 1 = 1.73 x 1(T 3 , 




-1000 -500 

V (MeV fm 3 ) 



FIG. 16. The dependence of 1Zr igSJ in SMEC on V in mir- 
ror systems 8 B(1^) (proton resonance) and 8 Li(l^) (bound 
state). 72-r using Eq. (|49|) (solid line) is compared to that us- 
ing expression (|50[) for the proton decay width (dashed line). 
The shaded area shows the range of IZr extracted in Ref. [l[ 
using experimental values of T p [52( and C v [53 [. 

whereas the experimental value extracted in Ref. |l| is 
(2.2 ± 0.2) x 10~ 3 . Results of SMEC obtained using 
Eq. ([50)) for the proton width become compatible with 
the experimental results for V < -1000 MeV fm 3 . The 
GSM prediction 1Z T = 1.88 x 10~ 3 given in Table |IX] 
is fairly close to experiment and to SMEC values. The 
SMEC values obtained with Eq. (pt9")l show a strong de- 
pendence on Vb ■ This estimate of TZr is valid in the range 
of < T p < T^, where T^) is the s.p. width, which 
in this case is 64.9 keV and 62.6 keV for p 3 / 2 and pi/2, 
respectively. 

The pi/2 and p 3 / 2 neutron contributions to the ANC of 
8 Li 1 + are shown in Fig.[T7J One can see that pi/ 2 and p 3 / 2 
contributions have different dependence on Vq which, in 
turn, reduces variations in 72-r ■ In general, for a physical 
range of Vq in this mass region, —500 MeV fm 3 > Vq > 
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FIG. 17. Squared ANCs of the P1/2 (a) and P3/2 (b) neu- 
tron contributions to Lij + . Shaded area mark experimental 



uncertainties on these quantities [l[. 



— 1200 MeV fm 3 , SMEC agrees somewhat better with 
experiment than SM. 



IV. OUTLOOK 

In the first part of the paper, we discussed the basic 
properties of ANCs and SPANCs. We broadly classified 
the behavior of SPANCS for charged and neutral parti- 
cles in terms of the Sommcrfcld parameter r\ and the wave 
number k, respectively, as well as the orbital angular mo- 
mentum. The extreme regimes of SPANC for charge par- 
ticles can be characterized by the complex turning point 
zt of the outgoing Coulomb wave function. We also dis- 
cussed the near-threshold behavior of ANCs. 

Based on the argument using the charge symmetry of 
the nuclear force, a simple relation (|30[) between proton 
and neutron ANCs in mirror pairs has been proposed (if . 
The estimate IZo is very useful as it relates cross-sections 
of low-energy direct and resonance proton capture reac- 
tions, which are difficult or impossible to measure, with 
neutron ANCs obtained in reactions with stable beams. 
In the second part of this study, the link between mirror 
ANCs through relation (|3TJ|) has been verified in our CSM 
calculations for different physical situations of the cou- 
pling to the scattering continuum and for various many- 
body states. It has been found that the key factor in 
determination of ANCs and the mirror ratio 1Z, as well 
as the sensitivity of ANCs to continuum coupling, is the 
distribution of spectroscopic strength that is both modcl- 
and interaction-dependent. For example, relative differ- 
ences of 1Z in GSM and VMC Q can be as large as 30%, 
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and the continuum coupling can change 1Z by up to 50% 
in exceptional cases. Also, differences with respect to 
IZo can be non-negligible. In this sense, ANCs and their 
mirror ratios are interaction-dependent. 

It has been found that the quantitative effect of the 
continuum coupling on ANCs and their ratios is minor 
if the spectroscopic strength is either localized in a sin- 
gle SM state or broadly distributed. This property is 
independent on binding energies of mirror states. On the 
other hand, if the spectroscopic strength is concentrated 
in several SM states, their coupling via the continuum 
space may result in a significant rearrangement of the 
spectroscopic strength; hence, appreciable variations of 
ANCs with respect to SM predictions. This effect is par- 
ticularly strong for near-threshold states that align with 
the decay channel. Since these special cases can be a pri- 
ori identified in standard SM calculations of the spectro- 
scopic strength distribution, the qualitative effect of the 
continuum coupling on SM results for ANCs and their 
ratios can easily be assessed without resorting to sophis- 
ticated CSM calculations, which ultimately provide the 
quantitative answer. 



Finally, let us state that uncertainties due to the model 
dependence can be significantly reduced if the effective 
SM interaction is optimized to both spectroscopic and 
reaction observables within a unified SCM framework. 
Work along these lines is in progress. 
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